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Abstract 

Action at a distance in Newtonian physics is replaced by finite 
propagation speeds in classical physics, the physics defined by the 
field theories of Maxwell and Einstein. As a result, the differential 
equations of motion in Newtonian physics are replaced in classical 
physics by functional differential equations, where the delay associ- 
ated with the finite propagation speed (the speed of light) is taken 
into account. Newtonian equations of motion, with post-Newtonian 
corrections, are often used to approximate the functional differential 
equations of motion. Some mathematical issues related to the prob- 
lem of extracting the "correct" approximate Newtonian equations of 
motion are discussed. 

Keywords: Delay equation, functional equation, inertial manifold, singular 
perturbation, radiation reaction, post-Newtonian 
AMS Classification: 34K60, 83C10, 83C25 



1 Introduction 

This paper — an expanded version of a lecture presented in the special session 
on Applied Dynamical Systems of the CAIMS 2001 meeting in Victoria, 

T Department of Mathematics, University of Missouri, Columbia, MO 65211. 
Email: carmen@chicone.math.missouri.edu 



1 



BC — is an invitation to explore some of the mathematical issues related to 
the foundations of post-Newtonian physics. 

Let us recall that Newtonian forces (for example, the inverse square law 
for gravitation) imply "action at a distance." This absurd, but outstandingly 
successful, premise of Newtonian theory predicts that signals propagate in- 
stantaneously. In classical physics, relativity theory postulates that signals 
propagate with a velocity that does not exceed the velocity of light. Thus, 
the forces of Newtonian physics must be replaced by force laws that take into 
account the finite propagation speed of the classical fields — electromagnetic 
and gravitational — which determine the forces acting on a moving body. In 
turn, the ordinary (and partial) differential equations of Newtonian physics, 
which are derived from the second law of motion d(mv)/dt = F, must be 
replaced by corresponding functional differential equations where the force 
F is no longer a function of just position, time, and velocity; rather, the 
classical force law must take into account the time-delays due to the finite 
propagation speed of the classical fields. 

The functional differential equations of motion for classical field theory 
are generally difficult, often impossible, to express in a form that is amenable 
to analysis. Thus, to obtain useful dynamical predictions from realistic mod- 
els, it is natural to replace the functional differential equations of motion by 
approximations that are ordinary (or partial) differential equations (cf. [^]). 
Of course, Newton's equations are the premier choice for approximating the 
true equations of motion. Indeed, due to the overwhelming success of New- 
tonian models in applied mathematics — where in most cases characteristic 
velocities are so low that relativistic effects are negligible — the dynamics of 
the true equations of motion are often ignored. The purpose of this paper 
is to discuss some of the mathematical issues that must be addressed to ob- 
tain a rigorous foundation for post-Newtonian dynamics, that is, Newtonian 
dynamics with relativistic corrections taken into account. 

I thank Bahram Mashhoon for suggesting several improvements of the 
original draft of this paper. 

2 Newtonian versus relativistic physics 
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2.1 Newtonian physics 



Let us consider the motion of a body p with mass m that is influenced by 
its Newtonian gravitational attraction to a second body P with mass M. In 
this case, the equation of motion for the body p is 

GMm (x — y) 
\x — yy \x — y\ 

where x denotes the position of p in space, y denotes the position of P, and 
G is the universal gravitational constant. Note that the force on the body 
p changes continuously with the position of P. By viewing changes in the 
gravitational force detected at p as a signal produced by manipulating the 
position of P, the (instantaneous) action at a distance of Newtonian gravity 
is seen to be equivalent to the infinite speed of propagation of gravity. 
By coupling equation flU) with the equation 

GMm (y-x) 

M y = —\ fa i r> 

\x — yy \x — y\ 

which models the motion of P as it is influenced by its gravitational at- 
traction to p, we obtain the prototypical model of Newtonian mechanics: a 
system of ordinary differential equations for the motion of two bodies influ- 
enced by their mutual gravitational attraction. The Newtonian model for 
two charged particles moving under the influence of the Coulomb force is 
essentially the same; the only change in the differential equations of motion 
is the replacement of the constant —GMm by the product of the charges of 
the particles. 

There are Newtonian models for systems of bodies, for fluids, elastic me- 
dia, etc. All of these models are ordinary (or partial) differential equations. 



2.2 Electrodynamics and gravitodynamics 

In classical physics, the fundamental electromagnetic and gravitational forces 
are generally given by time-dependent fields that propagate at the speed of 
light. 

Maxwell's field equations determine the properties of the electric and 
magnetic fields that influence the motion of charged particles through the 
Lorentz force law. In fact, the motion of a relativistic charged particle p 
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with mass m and velocity v, which is influenced by the electric field E and 
magnetic field B produced by a charged particle P, is given by 

d / mv \ . „ 1 

where c is the speed of light and q is the charge on the particle p. 

The electric field E that affects p at position x at time t is produced by 
P at y at the retarded time t — r where 

k(0 - y(* - r )l = cr 

(distance=speedxtime). Hence, the electric field at (x,t) is given by 

E(x,t) = f(t,x, T,y(t - r),y(t - r),y(t - r)), 

where / is some smooth function. The magnetic field has a similar form. 
It follows that the Lorentz force acting in spacetime at (x, t) is delayed by 
time r (itself an implicitly defined function of space and time) as p and P 
move, and therefore the equations of motion are a coupled system of retarded 
functional differential equations. 

Explicit representations of E and B are obtained from Maxwell's field 
equations 

div£ = 4vrp, cm\E = 

div B= 0, c curl B = 4nj + 

where p is the charge density, j is the current density, and c is the speed of 
light. In addition, we have the conservation of charge: divj = —dp/dt. 

A standard computation using vector calculus shows that the fields E 
and B are given by 

1 dA 

E = — — -grad0, B = cm\A, (3) 

where, once the Lorentz gauge condition 

divA+ = 
c at 

is imposed, the scalar potential <fi and the vector potential A satisfy the 
following wave equations with sources: 

1 d 2 ^ 



Id 2 A 4tt 
-AA= —j. 



dt 2 
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Using classical potential theory, the wave equation, and variation of param- 
eters, the retarded-time potentials are given by 

( )= fPjM^m dy , 

J \x-y\ 

A{ x ,t) = l - f ^-\^-y\/Q) d 

cj \x-y\ 

where the integrals are over all of space. The retarded-time electric and 
magnetic fields are then computed as in display (|3]). In practice, the current 
density is given by j = pv. Hence, by the conservation of charge, the charge 
density satisfies the continuity equation 

dp 

— + div(pv) = 0. 

The electric and magnetic fields produced by an arbitrary charge density 
are complicated. The potentials, however, can be computed explicitly for 
a point-charge moving in spacetime. In this ideal case, the charge density 
is a Dirac-type measure associated with the moving point-charge, and the 
resulting retarded-time potentials — called the Lienard-Wiechert potentials — 
are given by 

^ ' ^ ^ \x — y\ — y ■ (x — y)/c ^ re ^' 
A(x,t) = -<j)(x,t)y Tet , 

where the subscript ret indicates that y and y are evaluated at the retarded 
time t — t and r is given implicitly by the equation cr = \x(t) — y(t — r)\. The 
associated electric and magnetic fields are computed using the equations in 
display (^), and the corresponding equation of motion (|2]) for a point particle 
influenced by these fields is relatively simple. For example, the motion of 
two charged particles confined to move on a line, with only the Lorentz force 



taken into account, is modeled by the Driver- Travis system (see [10 , |18|) 

(1 - x2/c 2 )3/2 " C 2 miT 2 ( c _ (-iyxj(t - T tJ )) ' 
CTij \Xi Xj(t ^"ij)!' 

where is in the set {(1,2), (2,1)}. 
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There is a similar equation for the gravitational two-body problem (re- 
stricted to a line), but the equations are more complicated. The basic reason 
is that Einstein's field equation (R^ — {l/2)g llu R = StikT^, where R^ is 
the Ricci tensor, R is the scalar curvature, k = G/c 4 , is the stress-energy 
tensor, g^ v is the metric tensor, and the indices fi and v run over the integers 
from zero to four) is nonlinear. While the motion of a single test particle fol- 
lows a geodesic in spacetime, the general relativistic two-body system seems 
to be too difficult to write down explicitly. Models of this type are the sub- 
ject of current research. For example, relativistic effects, including radiation 
damping, are important in the dynamics of binary neutron stars (see and 
the references therein). In modern physics, the theoretical study of gravita- 
tional dynamics is generally more important than classical electrodynamics. 
The reason is that quantum mechanics has superseded classical electrody- 
namics, but there is yet no quantum theory of gravity. At any rate, the true 
equations of motion of the classical field theories are functional differential 
equations. 

3 Post-Newtonian approximation 

The basic idea of post-Newtonian approximation, from a mathematical point 
of view, is the expansion of model equations in powers of 1/c. From a physical 
point of view, the idea is to consider low velocity (compared with the speed 
of light) weak field limits. Note, for example, that the relativistic form of 
Newton's second law, where the rate of change of the momentum is given by 
|(mt;(l — \v\ 2 / (?)" 1 ' 2 ), reverts to Newton's law in the low-velocity limit. 

3.1 Radiation damping 

Classical mechanics is a mathematically consistent theory; it just 
doesn't agree with experience. It is interesting, though, that the 
classical theory of electromagnetism is an unsatisfactory theory 
all by itself. There are difficulties associated with the ideas of 
Maxwell's theory which are not solved by and not directly asso- 
ciated with quantum mechanics — R. Feynman [[12], p. 28-1]. 

According to Maxwell's field equations, a charged particle produces electro- 
magnetic fields as it moves. Since, in this particle radiates energy, 
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it must slow down. This basic intuition has led to some thorny issues in 
physics. To describe one of them briefly, let us consider a sphere consisting 
of identical charged particles. As this body accelerates, the various charges 
produce fields that affect the motion of the body through the Lorentz force. 

By considering motion along a line, Dirac reasoned that there are retarded 
and advanced self-forces that have (on average over the body) the post- 
Newtonian expansions 



ret 



aq z .. 2q 2 ... n . a . 
a& 3 c™ c 4 



„ aq 2 .. 2q 2 ... , a s 

where x is the position of the centroid of the sphere, a depends on the 
charge distribution (a = 2/3 for a round sphere), q is the charge on one of 
the particles, and a is the radius of the sphere (see fl2|, Ch. 28]). These forces 



are derived by expanding the Lienard-Wiechert potentials (see Ch. 9]). 

The existence of advanced forces (preaccelerations) would imply that the 
motion of a body is influenced by forces that are produced in its future. 
While this notion is surely problematic, it does lead to an interesting result. 

The ideal situation — Dirac's original motivation — is the theory of the 
electron, an elementary particle that seems to have no internal structure; 
it is supposed to be a point charge, that is, a sphere with radius zero. One 
manifestation of the internal inconsistencies mentioned in the quote by Feyn- 
man is the blowup of the coefficients of the first terms in these expansions 
as a — > 0. From Newton's second law, this coefficient, m e := aq 2 /(ac 2 ), 
represents a mass, called the electromagnetic mass. Thus, the difficulty can 
be stated as follows: the electromagnetic mass blows up as the radius of the 
sphere shrinks to zero. 

The radius of the sphere is given by a = a(q 2 /(m e c 2 )) and the number 
q 2 /{m e c 2 ) is called the classical electron radius (a is scaled out, because 
this coefficient depends on the shape of the original charge distribution). 
This value for the electron radius does not agree with the value that can be 
computed directly from relativity theory, a fact that led to a crisis in classical 



physics that is still not completely resolved (see [[12], Ch. 28] for an extended 
discussion). 

Dirac proposed a way to remove the apparent blowup of the electromag- 
netic mass. He theorized that the true self-force (on a moving electron) is 
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half the difference of the retarded and advanced forces; that is, 



F self = J™^( F adv- F ret)- 
Therefore, the self-force (the radiation reaction force) is 

2g 2 ... 
F self=3^- 

Using Dirac's theory, a post-Newtonian model for the motion of an elec- 
tron, confined to move on a line and with radiation reaction taken into ac- 
count, is given by the Abraham-Lorentz equation 

2g 2 ... „ 

mx = — - x +F, 

where F is some external force. In the presence of an electromagnetic field, 
the equation of motion in space would be 

1 x 2q 2 

mi ' . . > n 



q(E + -vxB) + -^i) + F. 



3c 3 

Since the particle radiates (produces fields that carry energy) the self force 
should cause the particle to lose energy and slow down. For this reason, 
the presence of the third-order time-derivative term in the first differential 
equation is called radiation damping. Is this intuition correct; that is, does 
the presence of this term cause damping? 

As a concrete example, consider two identical (isolated) charged particles 
each with unit mass that are one unit apart at rest. Imagine that these bodies 
are elastically connected by a force that can be modeled by Hooke's law with 
unit spring constant. Their relative motion, ignoring their mutual Coulomb 
and gravitational interactions, is modeled by the differential equation 

where e ~ c~ 3 . A similar scenario for two bodies with gravitational radiation 
damping taken into account, leads to a nonlinear oscillator of the form 



d 5 z 2 .. 
ez— — + z + z 
dt 5 
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where e ~ c -5 (see ||). In both cases, these post-Newtonian models are not 
Newtonian — for example, the differential equations are not second-order. 

Physical intuition suggests that the postulated two-body systems are 
damped oscillators. This leaves open the mathematical question: Is damped 
oscillatory motion predicted by these post-Newtonian models? 

The electrodynamic model is linear. Its characteristic equation 

-er 3 + r 2 + 1 = 

has the roots 

*-e + 0(e 2 ), -«-e + 0(e 2 ), l/e + 0(l); 

therefore, the first-order approximation of the general homogeneous solution, 
which is given by 

ae t/e + be~ et cos t + ce~ et sin t, 

has a "runaway" mode. A similar result is true, but more difficult to prove, 
for the nonlinear oscillator 

d 5 z 2 

ez—— + z + z = 1. 
dt 5 

Runaway solutions are clearly not physical. What do they represent? How 
should they be eliminated? More precisely, we may ask: What is the cor- 
rect Newtonian equation of motion with the radiation damping taken into 
account? 

There is an obvious mathematical answer to our question once our post- 
Newtonian models are recognized as singularly perturbed Newtonian equa- 
tions. By the definition of e as the reciprocal of a power of the speed of 
light, it is reasonable to assume that e is a small parameter, at least in the 
low-velocity regime. To recover the correct Newtonian model, we can apply 
Fenichel's geometric singular perturbation theory (see [JH], pj]). 

For the electrodynamic oscillator, we have the equivalent (singularly per- 
turbed) first-order system given by 

x = y, 

y = z, 

ez = z + x — 1. (4) 
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For the gravitodynamic oscillator, the appropriate first-order system is given 
by 



z — u, 

U = V , 

e^v = w, 
e 1/z w = x, 

e V3i = i-Z-V _ c l/3 5uX _ c 2/3 10vW 

2z 2 z z 

These systems are converted to regular perturbation problems by rescaling 
time. The corresponding fast-time electrodynamic first-order system is ob- 
tained by the change of variables s = t/e (s — t/e 1 / 3 for the gravitodynamic 
oscillator). These fast-time systems are equivalent to their original slow-time 
counterparts for e 7^ 0. In effect, the rescaling of time produces a new family 
of systems, still parametrized by e, but with a different limit as e — > 0. 
The fast-time electrodynamic system is given by 

x' = ey, 
y' = ez, 

z' = z + x — 1; (5) 

it is a regularly perturbed first-order system. The corresponding unper- 
turbed system (e = 0) has a two-dimensional invariant manifold, {(x,y,z) : 
z + x — 1 = 0}, consisting entirely of rest points. Moreover, this manifold 
is normally hyperbolic. In our special case, where the invariant manifold 
consists entirely of rest points, the requirement for normal hyperbolicity is 
that nearby trajectories are attracted to (or repelled from) the vicinity of 
the invariant manifold exponentially fast; or, in other words, the system ma- 
trix of the linearized system at each rest point has a zero eigenvalue with a 
two-dimensional spectral subspace (corresponding to the tangent directions 
on the invariant manifold) and a nonzero eigenvalue with a one- dimensional 
spectral subspace (corresponding to the normal direction to the invariant 
manifold). For the unperturbed system (|5|), the nonzero eigenvalue A is the 
same at each rest point. In fact, A = 1 and all nearby solutions are repelled 
from the unperturbed invariant manifold at this exponential rate. 

A fundamental result of Fenichel's theory states that a normally hyper- 
bolic invariant manifold persists. Hence, for each sufficiently small e 7^ 0, the 
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corresponding system (^) has a two-dimensional normally hyperbolic invari- 
ant manifold (not necessarily a linear subspace), called the slow- manifold. 
Moreover, the corresponding family of slow-manifolds depends smoothly on 
the parameter e and each slow-manifold is invariant, normally hyperbolic, 
and repelling. 

The original singularly perturbed slow-time family (Q) has a correspond- 
ing family of normally hyperbolic invariant manifolds. Each of them repels 
nearby solutions at an exponential rate. This accounts for the existence of 
runaway solutions. 

Using the invariance of the perturbed slow-manifolds, it is easy to ob- 
tain (in local coordinates) the equations of motion restricted to these two- 
dimensional manifolds. The perturbed invariant manifold, for the fast-time 
system, is the graph of a function of the form 

z = h(x, y) = 1 — x + ehi(x, y) + e 2 h 2 (x, y) + 0(e 3 ), 

where the hi are functions to be determined. Moreover, the dynamical system 
on this invariant manifold is given by 

x' = ey, y = ez = e(l — x + ehi(x, y) + 0(e 2 )). 

For e 7^ 0, the dynamical system on the corresponding slow-manifold for the 
original system is thus given by 

x — y, y = 1 — x + ehi(x, y) + 0(e 2 ). 

To solve for the functions hi, we use the invariance of the perturbed 
manifold. In effect, at each point on the manifold, each tangent vector X = 
(ey,ez,z + x — 1) corresponding to the system of differential equations (BI), 
is a linear combination of the basis vectors b x = (1, 0, h x (x, y)) and b y = 
(0, 1, h y (x, y)) for the corresponding tangent space on this manifold. The 
functions hi are obtained by equating coefficients of powers of e in the vector 
identity 

X = ab x + j3b y , 

where a and (3 are scalar variables (which are determined by the first two 
components of this identity). To first-order in e and in the original slow- 
time, the equation of motion on the invariant manifold is equivalent to the 
(Newtonian) second-order system 

i + ex + x = 1 , 
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a dynamical equation for a damped harmonic oscillator — just as it should 
be. 

For the gravitodynamic two-body system, a similar analysis results in the 
Newtonian dynamical equation 

15 

z + 32ezz(z — — ) + z = 1. 

Again, this is a (nonlinear) under-damped oscillator, at least if the initial 
separation between the bodies is near z — 1. 

3.2 Synthesis 

We have just seen that geometric singular perturbation theory (in particular, 
reduction to the slow-manifold) produces Newtonian model equations with 
post-Newtonian corrections that give physically reasonable dynamics; in par- 
ticular, the runaway solutions are eliminated. How can we justify using these 
models? Note, for instance, that the slow-manifolds in our models are unsta- 
ble; nearby solutions run away. In applied mathematics, we usually justify 
approximations by their stability. To validate the slow-manifold reductions, 
we must show that the resulting Newtonian model equations are "stable" 
with respect to the dynamics of the original functional differential equations, 
the true equations of motion in classical physics. 

Conjecture 3.1. In the low-velocity regime, a functional differential equa- 
tion of motion derived using the forces in classical field theory ( the Lorentz 
force or the relativistic gravitational force) has an inertial manifold X (a 
finite-dimensional, invariant, exponentially attracting, and smooth manifold) 
such that the restriction of the motion to this manifold is a Newtonian dy- 
namical system. The Nth-order singular perturbation problem, obtained by 
(post-Newtonian) expansion in powers o/l/c and truncation at order N , has 
an equivalent first-order system with a normally hyperbolic slow-manifold S. 
The corresponding vector fields that generate the dynamical systems onT and 
S agree to order N — 1 in 1/c. 

While the post-Newtonian expansion (and truncation) results in a system 
that might (and usually does) have runaway modes, these are simply arti- 
facts of the singular nature of the expansion. The long-term dynamics of the 
functional differential equations obtained by a direct application of classical 
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field theory is given by a Newtonian equation on a finite-dimensional inertial 
manifold, and the same dynamical system is obtained by appropriate reduc- 
tion to a slow-manifold for the singular high-order post-Newtonian equation. 
Runaway modes generally have no physical significance; equivalently, the 
slow-manifold is generally not an attractor. 

The conjecture states a rigorous justification for the validity of the post- 
Newtonian model that agrees with the dynamics on an inertial manifold 
which is supposed to exist in the low-velocity regime. But, what is a low 
velocity? A deeper question addresses this issue directly. For what range of 
parameter values does the inertial manifold persist? 

In the mathematical analysis, 1/c, or (even better) a characteristic veloc- 
ity divided by c, is to be viewed as a small parameter. If the conjecture is 
valid, then there is a lower bound for e such that the corresponding functional 
differential equation has an inertial manifold. A theorem meant to validate 
the conjecture should include an estimate for this bound. 

What accounts for the transition from the low velocity to the high velocity 
regime? Answer: As a characteristic velocity increases, a bifurcation will 
occur that destroys an inertial manifold and, therefore, the validity of post- 
Newtonian approximation. How can such bifurcations be detected? 



4 Delay equations 

To test the conjecture stated in the last section, let us replace the functional 
differential equations of mathematical physics (with space-dependent delays) 
with families of delay differential equations of the form 

±(t) = f(x(t),x(t-r)), x G R n , tgR, (6) 

where the delay r replaces the small parameter corresponding to 1/c and / 
is a smooth function. 

The usual state space for the delay equation (Q) is C([— r, 0]), the space 
of continuous functions that map the interval [— r, 0] into R n . A basic result 
in the well-developed theory of delay equations (see || |9|, [El], states that 



the initial value problem consisting of equation (|6]) and an initial function in 
C([— r, 0]) has a unique solution t \— > x(t) for t in some interval [0, /?), where 
(3 > or (3 = oo. For such a solution, the state of the system at time t > 
is given by x t G C([— r, 0]), where x t (9) = x(t + 6). This assignment defines 
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a semi-flow in the state space given by (T*-?/>)(#) = x(t + 9), where t h- > 
is the solution with initial condition if>. 



4.1 Inert ial manifold reduction 

A heuristic argument indicates that the delay equation has an inertial 
manifold for small |r|. Let us first introduce a new variable t = st so that 
y(s) = x(st) is a solution of the delay equation 

y(s) = rf(y(s),y(s-l)) (7) 

whenever x is a solution of the delay equation (|6]). The state space for the 
family (|7|) is C([— 1,0]). Also, the unperturbed system (r = 0) 

y(s) = 

generates the semi-flow given by 
TV 

The n-dimensional submanifold 



(ip(t + 9), < t < 1 and t + < 0, 
"^(O), otherwise. 



J : = e C([-l, 0]) : ^(0) = a for some a G R n } 

consists entirely of rest points and is normally hyperbolic. In fact, every 
solution reaches X in time £ = 1, a rate of normal contraction that is faster 
than any exponential decay. Thus, it is reasonable to expect that X persists 
as an inertial manifold for sufficiently small \t\. 

We have not proved the persistence of the invariant manifold. Because the 
state space is infinite-dimensional, there are some delicate issues involved in 
proving the existence and smoothness of an inertial manifold. Recent results 
(see, for example, |l]) on the persistence of infinite-dimensional normally 
hyperbolic invariant manifolds may apply. But, at the first level, a theorem in 
this direction would state the existence of an inertial manifold for sufficiently 
small values of the parameter. This would not be completely satisfactory for 
future applications to physics where some estimate of the relevant parameter 
values would be required. The foundation for a direct proof of the existence 
of an inertial manifold, with explicit bounds but without a proof of the 
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required smoothness, is contained in the work of Yu. A. Ryabov and R. D. 
Driver (see [7|, ||). The smoothness of the inertial manifold is proved in ||. 

Having indicated that an inertial manifold exists, let us reduce the dy- 
namical system to the inertial manifold and thus obtain the "Newtonian 
system" corresponding to the delay equation. To do this, suppose that £ is 
a coordinate on R n and y(t, £, r) is the flow on the inertial manifold; that is; 
t i — ► y(t,£,r) is the solution of the delay equation such that y(0, £, r) = £. 
More precisely, ip{9) = £ is the initial condition for the solution y of the delay 
equation. With this notation, we have that 

y{t,M =f(y{t,Z,T),y{t-T,Z,T)); 

hence, the vector field that generates the flow Y on the inertial manifold is 
given by 

X(£,r):=y(0,e,r) = /(£,y(-r,£,r)). 
Its expansion at r = is 

X(£, r) = /(£, - rD 2 f(t + ^ £)(/(£, 0, fit, 0) 

+ 0(A/(£, + 3£> 2 /(e, 0)/(e, 0)} + 0(r 3 ), (8) 

where the operator respectively D 2 , denotes differentiation with respect 
to the first, respectively the second, argument of /. 

It is interesting to note that the presence of a small delay in a conservative 
system often results in damped long-term dynamics on an associated inertial 
manifold. For example, the Duffing-type model equation 

x + tu 2 x = —ax(t — r) + bx 3 (t — r) 

with small delay r in the restoring force, reduces (by a formal computation 
to first-order in r) to the van der Pol-type model equation 

x + r(3bx 2 — a)x + (a + uj 2 )x — bx 3 = 

on its inertial manifold. This example illustrates a phenomenon that is remi- 
niscent of quantization: while most periodic solutions in one-parameter fam- 
ilies of periodic solutions in a conservative system disappear in the presence 
of a small delay, some persist as limit cycles. Does this observation have 
physical significance? 
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4.2 The analog of post-Newtonian expansion 

For the delay equation (|6]), the analog of post-Newtonian expansion is the 
expansion of the function r i— > f(x(t),x(t — r)) at r = 0, where the first few 
terms of the series are given by 

f(x(t),x(t - r)) = f(x(t),x(t)) - rD 2 f(x(t),x(t))x(t) 

+ T -(D 2 f(x{t),x{t))x{t) + D 2 J(x(t), x(t))(x(t), x(t))) 
+ 0(r 3 ). 

Truncation of the expansion at order N in r produces an iVth order ordinary 
differential equation of the form 

(-lf^D 2 f(x,x)x^ =F{x,x,... (9) 

the desired analog of a post-Newtonian expansion in classical physics. More- 
over, by setting /i := t 1 ^ n ' 1 ^ and treating /i as a small parameter, we obtain 
a singularly perturbed first-order system of the form 

N • 

V 2/i = 2/2, 
H N y N -2 = Vn-i, 

fM N y N -i = F(x, yi, ... , y N -i, /J, N ~ 1 ). 

Under the assumption that D 2 f(x, x) is invertible, the geometric singular 
perturbation theory can be applied to prove that this system has an n- 
dimensional slow-manifold for sufficiently small r^O. 

4.3 Synthesis 

For delay equations, the family of flows restricted to the family of slow- 
manifolds (parametrized by r) is generated by the family of vector fields 
F(£,r). A result in H] states that this family of vector fields agrees with 
X(£,t), the family that generates the flows on the inertial manifolds, to 
order r 2 . 
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For applications to physics, the result that X and Y agree to order- 
two in t is sufficient for many applications. Of course, it is not difficult to 
show that these vector fields agree to order three, four, etc. On the other 
hand, the complexity of the terms that appear in the expansion increases 
rapidly with the order. While the equality of the low-order coefficients of 
t in the two expansions can be proved by direct computation, there does 
not seem to be an easy abstract proof for the equality of all coefficients. 
Thus, the conjecture that X and Y agree to order N — 1 (one order less 
than the order of the "post-Newtonian" truncation) is open. The hypothesis 
that D 2 f(x,x) is invertible should be sufficient for the conjecture, at least 
in the case where N is sufficiently large. Under the stronger hypothesis that 
D 2 f(x, x) is infinitesimally hyperbolic (no eigenvalues on the imaginary axis), 
no such restriction on N should be necessary. 

A apparent difficulty to be overcome in the proof of the conjecture for 
delay equations is encountered in the proof of the conjecture for the special 
case of linear delay equations (see ||). 

Theorem 4.1. If A is invertible, |r|||A||e < 1, and 

x(t) = Ax(t - r), 

then the family of vector fields that generates the flow of this system on its 
family of inertial manifolds (parametrized by r) agrees to order N — 1 with 
the family of vector fields on the corresponding family of slow-manifolds for 
the corresponding family of Nth order systems (j^). Moreover, the family of 
vector fields on the inertial manifolds is given by 



This nontrivial identity can be proved using Abel's generalization of the 
binomial theorem; namely, the identity 




The proof in || requires the combinatorial identity 




a 




(a + (3)(a + (3 + m) 



m—l 
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(see [|T3"| , p. 19]). It seems that a "nonlinear" replacement for this combina- 
torial identity will be required to prove the analog of Theorem |4.1| for the 
delay equation (^). 
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